Self-organized metal nanostructures through laser driven thermocapillary convection 
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When ultrathin metal films are subjected to multiple cycles of rapid melting and resolidification 
by a ns pulsed laser, spatially correlated interfacial nanostructures can result from a competition 
among several possible thin film self-organizing processes. Here we investigate self-organization 
and the ensuing length scales when Co films (1 — 8 nm thick) on SiC<2 surfaces are repeatedly and 
rapidly melted by non-uniform (interference) laser irradiation. Pattern evolution produces nanowires, 

which eventually break-up into nanoparticles exhibiting spatial order in the nearest neighbor spac- 
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ing, \nn2- For films of thickness ho > 2nm, Xnn2 <x h ' while the particle radius varies as 

1/2 

r P 2 oc h ". This scaling behavior is consistent with pattern formation by thermocapillary flow and 
a Rayleigh-like instability. For h < 2 nm, a hydrodynamic instability of a spinodally unstable film 
leads to the formation of nanoparticles. 

Metal nanostructures, possessing spatial order and size uniformity are attractive for a variety of ap- 
plications including: optoelectronics, plasmonics, chemical and biomedical sensing, and catalytic devices 
IlL l2l La, |4j, |5Q. In recent years there has been tremendous research emphasis on finding self-organization 
strategies to make strongly correlated metal and semiconductor nanostructures. In the area of thin film 
pattern formation, epitaxial strain driven self-organization in crystallographic systems has been somewhat 
successful 1 6, 7, |1 9]. However, this strategy of exploiting epitaxial strain cannot be applied to materials on 
important amorphous surfaces, such as Si02- One possibility for self-organization on amorphous surfaces 
is via dewetting instabilities, such as those observed in liquid polymer films on various inert surfaces, in 
which preferred length scales appear when destabilizing long range attractive intermolecular forces over- 



come stabilizing surface tension forces 



12]. We have previously shown that multiple pulses of 



single beam ns laser irradiation of Co films with thickness 1 < ho < 8 nm on Si02 surfaces demonstrate 
self-organization by a thin film hydrodynamic (TFH) instability with length scales characterized by the 
metal surface tension (7), the dispersion force characterized by the Hamaker constant A and film thickness 
ho Il3ri . The ensuing pattern had short range spatial order (SRO) in the nearest-neighbor (NN) particle 
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spacing Xnni, which scaled with thickness as h^, while the particle radius varied as r p \ ~ h^ 3 . R ecent ly 
we have also provided preliminary evidence that 2-beam laser interference irradiation can give rise to quasi- 
2D ordering of the nanoparticles with both long range order (LRO) and SRO MM . However the evolution 
of the 2-beam pattern from the flat film to nanoparticle state was not investigated, and the mechanisms for 
spatial ordering were not identified. Here we report detailed experiments of 2-beam interference irradiation 
of Co on Si02 as a function of film thickness for 1 < ho < 8 nm that demonstrate that the pattern evolution 
and resulting length scales are consequences of different mechanisms, including the TFH instability, ther- 
mocapillary (TC) flow and a Rayleigh-like instability, based on the timescales of these processes. One can 
access a particular mechanism via the choice of film thickness and/or thermophysical parameters such as 7, 
A and the temperature gradient. 

Cobalt films ranging in thickness from ~ 1 — 8 nm were deposited by e-beam evaporation onto optically 



ito optic 
lill5 



smooth Si02 surfaces in vacuum. Specific details of the experiments are published elsewhere 1141 11511. In 
brief, following deposition, the Co films were irradiated in vacuum with a 266 nm pulsed laser operating 
at a repetition rate of 50 Hz with a pulse time r ~ 9 ns by single beam I13I or 2-beam interference 
irradiation. In the single beam case, films were exposed to an unfocused beam at normal incidence. For the 
2-beam condition, an interference angle of ~ 45° generated a periodic laser intensity profile whose contrast 
was maximized by adjusting the energy density of the off-normal beam with a 750 mm lens positioned 
300 mm away from the sample at angle of ~ 45°. In both instances, Co films were irradiated for similar 
times as measured by the number of laser pulses, n, for 5 < n < 10, 500. Also, comparable laser energy 
densities were used for both types of irradiation, with the energy density, Ei aser between Exhreshoid < 
60 < Ei aser < 200 mJ llql . The resultant morphological evolution was investigated as a function of n and 
was examined and characterized using a Hitachi S-4500 scanning electron microscope (SEM). It was also 
verified that all the films studied had minimal evaporation, as confirmed by performing energy dispersive 
X-ray spectrometry measurements of the Co concentration after irradiation. 

As detailed in prior studies for single beam irradiation HI 16], the film shows several distinct patterns 
enroute to a final robust state of nanoparticles characterized by a spatially correlated nearest neighbor (AW) 
spacing. Here we briefly summarize the trends for the purpose of comparison with the pattern evolution 
under 2-beam interference irradiation. More detailed analysis the single beam data can be found elsewhere 
ll3,Q]. At the early stages (n = 10) regular holes with a characteristic diameter form, followed by a 
cellular network of polygonal structures. Eventually, the polygonal networks evolve into nanoparticles, 
which form predominantly at the vertices of the polygons. An important property of the evolution is that 
the patterns possessed a characteristic length scale at every stage. In Fig. [Ha), the scaling behavior of the 
spatial correlation in the final stable nanoparticle state is shown for varying initial film thicknesses. The 



observed trend with ho was in agreement with classical linear TFH dewetting theory 



11] in which: 



A TFf , = V ±^ (1) 

where A-tfh represents the average NN nanoparticle spacing, Xnni, A = 1.4 x 1CP 18 J is the experi- 
mentally estimated Hamaker constant for the SiO-y] Co /Vacuum system , 7 = 1.88 ^ is the Co liquid 
surface tension and ho is the film thickness | u], X1TI . Furthermore, volume conservation implied that the 



nanoparticles will have a radius that will vary with the initial film thickness ho as: 
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where f(9) accounts for the contact angle of the nanoparticles 11311 . We experimentally observed a 
monomodal particle size distribution with the average radius shown in Fig. (Hb) and a trend consistent 
with the above theoretical prediction from linear TFH theory. 

The results of 2-beam irradiation were qualitatively different from single beam irradiation. In Fig. [2ja- 
c), the typical pattern evolution for films with ho > 2 nm is shown as a function of n. In this case, the 
early stages are comprised of spatially periodic film rupture at length scales comparable to the interference 
spacing. Longer irradiation yielded the formation of long, cylindrical-like "nanowires" and continued irra- 
diation resulted in the break-up of these nanowires into particles. The final particle state is characterized by 
a quasi-2D, comprised of the LRO due to periodic laser intensity and SRO resulting from the break-up of 
the nanowires. The observed differences from the two processing conditions can be explained on the basis 
of the various mechanisms of fluid motion operating under the two irradiation conditions. When irradiating 
the metal film with a 2-beam interference pattern, the resulting periodic laser intensity induces a transient 
and periodic thermal gradient along the plane of the film which creates a surface tension gradient not present 
in the single beam irradiation. Consequently, TC or Marangoni convection of the molten Co can occur. To 
contrast the results of single beam and 2-beam irradiation, the time scales for the various mechanisms were 
estimated. The timescale for TC flow can be expressed as 
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where A; aser = 350 nm is the laser fringe spacing, n = 4.45 x 1CT 3 Pas, is the Co metal viscosity, 
t^, = —0.34 x 10~ 3 Jm~ 2 K~ 1 is the rate at which the surface tension of the Co changes with temperature 
and AT is the maximum temperature difference between the peak and valley of the laser fringe. The 
timescale for the TFH instability can be expressed as |Q, 11 ] : 
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An important aspect of the above equations is that the time scale of the TC flow decreases with increasing 
ho while the TFH instability increases with increasing ho- We determined that typical liquid lifetimes range 
from 1 < tl < 10 ns and thermal gradients range from 0.5- 0.7 J- flfl. From Eq. [3]and|4j t s > tmo, 
for films with h D < 2 tim and t s < tmo, f° r films with ho > 2 tim. This implies that pattern formation 
is dominated by the TC flow for ho > 2 nm with the TFH flow dominating for ho < 2 ran, as shown in 
the comparative timescale plot in Fig. [3j This result was confirmed experimentally as nanowire formation 
was observed only above 2 ran while it was absent for films < 2 nm, as shown in Fig. Oinset A). For film 
thicknesses near the cross-over point of ho = 2 nm, the pattern consisted of particle formation with some 
evidence for lateral movement of the metal (Fig. [3]lnset B)) indicating that both mechanisms are operative 
111. 

For films with ho > 2 nm nanowire formation under the 2-beam irradiation also permitted access to a 
Rayleigh-like instability in which cylinders are unstable to wavelengths > 2nr cy i, with the fastest growing 
wavelength scaling as A# oc r cy i , where r cy i is the radius of the cylinder. The characteristic time scale for 



this process can be expressed as u2 , 



l P r lyl 



(5) 



where p = 7.8 is the liquid density of Co [20]. Assuming these are partial cylinders with contact angles 



> 90°and using the projected width of these regions from SEM images, we estimated r cy i, preceding 
the break-up into nanoparticles, to be in the range of 35 — 55 nm. Based on these values of r cy i, the 
typical magnitude for the cylinder break-up time was estimated as TR ay ~ 1 — 2 ns. Since this time scale 
is comparable to or smaller than the typical liquid lifetimes, the Rayleigh process is clearly accessible 
whenever a cylinder is formed. Also, according to the Rayleigh formulation, the particle radius scales 



21, 



22J]. Through 



linearly with both the spacing between particles and the radius of the original cylinder 11191 
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volume conservation, the cylinder radius varies as r cy \ oc h , resulting in a similar scaling relation for 
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the size of the particle as r P 2 oc hj ■ As shown in Fig. QJa-b), the particle spacing and the size both 
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scale as hj and the ratio of particle spacing to radius is independent of film thickness. While, the ratio 



in the classical Rayleigh break-up of a perfect cylinder is expected to be 



AjVJV2 
r p 2 



4.7 (12, 



21, 



22] our 



measured ratio was ~ 5.6. We have ruled out previously existing explanations for how this ratio can be 



modified [23, 



24, 



25, 



26, 



27, 



28|]. Nichols et al. [23] have shown that if the break-up is dominated by solid- 
state surface diffusion, then the ratio is maintained, however, if the flow is generated by external volume 
diffusion, then the ratio grows to ~ 6.1. However, the timescale over which external volume diffusion 
can cause break-up has been estimated as ~ 0.3 s (t vo i ~ 13 fkT^ ' wh. ere D v is the volume self-diffusion 



coefficient, Q is the atomic volume and / is the correlation factor 12311 ) for cylinder with r cy i = 40 nm. 
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This timescale is well outside the realm of the laser processing time (for t = 0.3 s, n ~ O(10 6 )), and 
as such can be disregarded as a possible break-up mechanism. We have also estimated the correction to 
the ratio for situations in which truncated cylinders break-up based on the work by McCallum et al. 
and determined that only a small change to the ratio will be introduced. While there have also been papers 
examining the effect of having a non-uniform contact line or an anisotropic surface, these models also do 



not yield the observed 5.6 ratio 127 



28]. It is possible that the thermocapillary flow and a Rayleigh-like 
instability can drive self-organization of the film such that the liquid profiles are stationary solutions to the 
non-linear equations describing the pattern evolution |2^, 3(1. However, a nonlinear analysis along with 
further experiments are required to verify such a hypothesis. 

In conclusion, self-organization leading to spatially correlated nanostructures under ns laser irradiation 
of ultrathin Co films on Si02 has been investigated. For films with thickness 2 — 8 nm, non-uniform laser 
irradiation by 2-beam interference leads to pattern formation characterized by the formation of nanowires 
via thermocapillary flow and eventually to nanoparticles via a Rayleigh-like break-up of the nanowires. 

1/2 

The nanowire break-up results in an average nearest-neighbor particle spacing of Xnn2 °c h with radius 
r p 2 oc h^/ 2 in contrast to the TFH instability in which Xnni °c and r p \ oc h^ 3 . For films with 
ho < 2 nm nanowire formation was absent because the TFH time scales were much shorter than the TC 
timescale. These results show that self-organization through laser-induced hydrodynamic flow can be used 
to make a variety of strongly correlated surface nanostructures on amorphous substrates by selecting the 
appropriate film thickness based upon thermophysical parameters. 
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Figure captions 

• Figured] Self-organized length scales versus initial film thickness: (a) \/Kq dependency for the parti- 
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cle radius from 2-beam interference irradiation (open squares); and dependency from single beam 
irradiation (closed triangles), (b) \fho dependency for the particle spacing from 2-beam irradiation 
(open squares) and dependency from single beam irradiation (closed triangles). 

• Figure [2j SEM images depicting the stages of pattern formation for 2-beam interference irradiation 
of a ~ 6 nm film as a function of increasing number of laser pulses n: (a) periodic rupture (b) 
nanowires and (c) the final nanoparticle state which exhibits both LRO and SRO. (Co rich and Si02 
rich regions correspond to bright and dark contrast respectively). 

• Figure [3j Comparison of the time scales for TFH dewetting and thermocapillary (TC) flow as a 
function of initial film thickness with the morphology for various film thickness following irradiation 
(insets). (A) h < lnm film showing that the entire film dewets; (B) h ~ 2 ram film showing that 
while TC flow causes the film to split there is still evidence for TFH dewetting and (C) h ~ 4.5 ram 
film showing that TC forces completely dominate morphology change. (Co rich and Si02 rich 
regions correspond to bright and dark contrast respectively). 
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